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SOME MEASUREMENTS OF ATMOSPHERIC TURBULENCE 
OBTAINED FROM FLOW-DIRECTION VANES 
MOUNTED ON AN AIRPLANE 


By Robert G. Chilton 
SUMMARY 


The power spectrum of relatively short wavelength turbulence in 
the atmosphere was calculated from measurements made in flight. The 
range of wavelengths covered by these measurements was from 10 feet to 
200 feet. The power spectral density varies with the square of the 
wavelength of the turbulence. This variation is in agreement with the 
high-frequency asymptote of the spectrum form generally assumed for 
isotropic turbulence. Flow-direction vanes were used to measure the 
vertical and horizontal components of gust velocity normal to the 
flight direction. The power spectral densities of the two components 
are, for practical purposes, equal. The use of vanes is shown to 
afford a simple, direct method of obtaining the power spectral density 
of atmospheric turbulence in the relatively high frequency range of 
airplane response. 


INTRODUCTION 


The introduction of the methods of generalized harmonic analysis 
to measurements of atmospheric turbulence has provided a form of sta- 
tistical information which makes it possible to analyze problems con- 
cerning the motion or the flight path of an airplane which is flying 
in rough air. Because of the random nature of turbulence, previous 
methods available to airplane designers for describing rough air pro- 
vided only statistical estimates of peak values of gust velocity, 
velocity gradients, and the frequency of their occurrence. Although 
this information is useful for calculating maximum values of airplane 
gust response, it does not describe the gust disturbance in sufficient 
detail for calculating the average effect of the continuous gust dis- 
turbance on the response of an airplane which is encountering atmospheric ks 
turbulence. The &irplane response to gusts is determined by the dynamic 
characteristics of the airplane. The response is, therefore, dependent 
upon the gust wavelength. Calculating the response of an airplane in 
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turbulent air requires information pertaining to the gust disturbance 
which distinguishes between gusts of different wavelengths. The appli- 
cation of the methods of generalized harmonic analysis to gust data 
accomplishes this purpose. 


The function which gives the information іп the desired form is 
called the power spectral density. The power spectral density of gust 
vertical velocity, for example, is & function of the inverse wavelength 
which gives a measure of the mean-square velocity of all those compo- 
nents of vertical velocity having wavelengths within an incremental 
band width. The mean-square gust velocity is, therefore, equal to the 
area under the power-spectral-density curve. The mean-square value of 
a function has by electrical analogy come to be considered in terms of 
average power. The power spectral density of a function, therefore, 
describes the manner in which the total average power of the function | 
is distributed over the entire frequency range. A complete mathematical 
development of the concepts of the power spectral density of a random 
function can be found in a number of texts (for example, ref. 1). Some 
particular applications of the power-spectrum concept to airplane 
response to atmospheric turbulence can be found in references 2 and 3. 
some of the problems involved in the calculation of the power spectrum 
from a finite time series are treated in reference Ш. 


The usefulness of the power-spectrum concept of random functions 
lies in the theorem that the power spectrum of the response of a linear 
system subjected to a random disturbance is equal to the power spectrum 
of the disturbance multiplied by the amplitude squared of the frequency- 
response function of the system. It is this relationship which makes 
possible the evaluation of the effect of atmospheric turbulence on air- 
planes and airplane control systems. Some examples of problems con- 
cerning atmospheric turbulence which can be analyzed in this manner 
are airplane passenger comfort, tracking or bombing errors due to flight- 
path response of airplanes, and the effects of the combined spanwise 
and chordwise variations of angle of attack on a wing in turbulent air. 


Before these methods can come into general usage, the designer 
mist have available reliable information in the form of spectral densi- 
ties of atmospheric turbulence. The need exists, therefore, for exten- 
sive measurements of turbulence over a wide range of wavelengths repre- 
senting many different meteorological and topographical conditions. 
With sufficient data, it can be determined to what extent the nature of 
turbulent air can be generalized or what special conditions must be 
imposed before a generalized power spectrum can be applied to specific 
problems. The first work in this direction (ref. 2) presents data and 
describes & method of obtaining the power spectrum of gust velocities 
in the frequency range normally called the short-period range of an  . 
airplane. For the airspeed of the airplane in reference 2, this short- 
period range corresponds to & range of gust wavelengths from approximately 
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200 feet to 2,000 feet. Other work has been done to extend measurements 
to very low frequencies or long wavelengths. 


The purpose of this paper is to present the spectrum of some short- 
wavelength measurements of gust velocities. The range of wavelengths 
covered is approximately from 200 feet to 10 feet. The data presented, 
when considered with measurements at the lower frequencies, are useful 
in extending to higher frequencies the data available on which approxi- 
mations of gust power spectral densities can be based. The measurements 
were obtained by the use of flow-direction vanes mounted on & boom 
ahead of the nose of an airplane. This method provides a ready means 
for accumulating a large amount of data in this frequency range and it 
can be employed simultaneously with other methods for making measure- 
ments in the lower frequency ranges. 


SYMBOLS 
G(r) power spectral density 
f frequency, cycles per second or cycles per foot 
T finite interval of time, sec; also, scale of turbulence, ft 
w vertical velocity, ft/sec 
t time, sec 
1 distance from angle-of-attack vane to center of gravity of 
airplane, ft 
j = -i 
X longitudinal axis of reference fixed in airplane 
2 normal axis of reference fixed in airplane 
8 airplane angle of pitch, radians 
V airplane forward velocity, ft/sec 
a angle of attack, radians 


[ 1 frequency-response-function notation 
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Subscripts: | ; & 2 = 
ай input 

о output 

É gust 

v vane-indicated - 
a airplane 


A dot over a quantity indicates differentiation with respect to time. 
INSTRUMENTATION AND TESTS 


Flow-direction vanes on & fighter airplane were used to record 
angle of attack and angle of sideslip (fig. 1). The vanes were mass 
balanced and located 54 inches ahead of the nose of the airplane. The 
undamped natural frequency of the vanes at sea level and at a test 
airspeed of 220 feet per second was approximately 45 cycles per second. 
The damping ratio was 60 percent of critical. These values were 
obtained by wind-tunnel tests. The recording system was properly 
damped and had an undamped natural frequency of approximately 100 cycles 
per second. Standard NACA airspeed-aitimeter, linear-acceleration, 
angular-velocity, and attitude-pitch recorders were used. A statoscope 
having a sensitive differential pressure cell with a scale range of 
£5 inches of water was employed to measure variations in pressure 
altitude. 


The airplane was flown in a straight path at an altitude of approxi- 
mately 1,000 feet with a minimum of control application for periods of 
1 minute or less. Test runs were made with the airplane flying both 
parallel and perpendicular to the wind direction over commercial and 
light-industrial areas in the vicinity of Portsmouth, Va. The test 
runs were made in clear air on two afternoons in late summer. The 
dates and other information pertinent to the test runs are given in 
table I. 


METHOD OF ANALYSIS 


Definitions 


The relationships which pertain to the concept of the power spectral 
density of a random function are presented here as defined in reference 1l. 
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Qust vertical velocity is chosen as the function for which the spectral 
density is of interest. If ча(%) is defined as the vertical component 
of gust velocity, the spectral density of Wg can be defined in the 


following manner. For -T St < Т, let 


f(t) = wg(t) (la) 
and elsewhere, 
f(t) = O (15) 


The Fourier transform of f(t) exists and is by definition 
оо T 
F(f) -/ £(t)e" dt tat -/ £(t)e" 250 %а+ (2) 
=O ~ 
The power spectral density of gust vertical velocity is then defined as 


2 
G(f) = a = (е) (3) 


for which only positive frequencies are considered. The mean-square 
value of gust vertical velocity Wee is found by integrating the power 


spectrum 
2 
Wg - а(ғ)ағ (4) 
O 


The output spectrum Со(Ғ) of the response of an airplane (which is 
assumed to be a linear system) in rough air which has the spectrum G4 (f) 
is (ref. 1, p. 288) 


Golf) = |Ү(2х3ғ) |204 (£) (5) 


where Y(2rjf) is the frequency-response function which relates the 
response of the airplane to the gust disturbance. 
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Measuring Technique 


Most of the methods of obtaining the power spectrum of atmospheric | > 
turbulenée which have been reported require the use of.the frequency- 
response function which relates the response of the airplane to the f 
gust input. In utilizing this technique, the response of the airplane | m 
in rough air is measured and the power spectral density of the response 
is computed. The power spectrum of the gust-velocity input is then 
calculated by using equation (5). When possible, the response is 
measured with the controls fixed. If the length of record required is 
too great for the airplane to stay in trim without control corrections, 
the effects of the control inputs must be removed and the analysis 
becomes cumbersome. In this event, a stabilized airplane may be used, 
in which case it is necessary to know the frequency-response function 
of the airplane-autopilot combination. E 


Since the measurements described in this paper are confined to 
high-frequency gusts (1 cycle per second and higher at an airspeed of 
220 feet- per second), it is possible to employ a flow-direction vane 
as the primary means of measurement. The gust velocity can be calculated 
from the indications of the vane and corrected for the effect of the 
airplane response to the gusts. No frequency-response functions for 
the airplane are required, and no restrictions on the application of 
controls are necessary. 


Geometrical Relationships 


The geometrical relationships for the airplane and flow-direction 
or angle-of-attack vane are shown in figure 2. All angles are assumed 
to be small variations from the steady state, and small angle values 
are assumed for the trigonometric functions. If variations in upwash 
are neglected, the vane-indicated angle of attack a, can be expressed 
as 


Iz 


Ол = Qa + ас - T (6) 


where 1 is the distance from the vane to the center of gravity of the 

airplane. The term 10 is the vertical velocity &t the vane resulting 

from pitching velocity of the airplane. The boom is assumed to be rigid. 
W Wo. 

Since a, = 9 + = and Qg = - the gust velocity Wg can be written Е 


as 


Wg = V(ay - 9) - Wa + 16 (7) 
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The only quantity not measured directly is the airplane vertical velocity 
which can be calculated by integrating the normal accelerations as dis- 
cussed in the following section. 


Method of Calculation 


The quantity (gust vertical velocity) for which the power spectral 
density is required is given in equation (7) as the sum of measured or 
calculated quantities. The primary quantity is the vane measurement 
and those describing airplane reactions are correction terms. The 
simplest way of handling this problem is.to perform first the arithmeti- 
cal calculations indicated in equation (Т) to obtain a single record of 
computed gust vertical velocities. The results can then be processed 
to obtain the power spectral density of the gusts. Тһе only quantity 
not measured directly is the vertical velocity of the airplane. Normal 
acceleration and variations in pressure altitude were measured to pro- 
vide information for calculating the airplane vertical velocity. The 
normal accelerations were integrated and the slope of the pressure- 
altitude curve provided the initial value of vertical velocity. The 
calculated vertical velocities were integrated a second time for com- 
parison with the pressure-altitude curve. 


Since gusts of relatively short wavelength were of primary interest, 
the assumption that the angle-of-attack vane gives a sufficiently accu- 
rate measurement of the gust vertical velocity without corrections 
(w; = Vay) was considered. This assumption would provide a great saving 
іп the film reading required to correct for the airplane response. For 
this purpose & trial record of approximately 20 seconds duration was 
analyzed to determine the importance of the corrections. The power 
spectral density of the gust vertical velocity calculated from equa- 
tion (7) is compared in figure 3 with that obtained by assuming no 
airplane response. The comparison indicates that, above 1 cycle per 
second, it is unnecessary to correct the vane measurements for the 
effects of the airplane response to the gusts for the particular airplane 
used. The vertical-velocity response of the test airplane in the range 
of frequencies of interest is so small that the corrections may be neg- 
lected. The frequency range for which this effect is true for other 
airplanes depends, of course, upon the dynamic characteristics of the 
particular airplane. 


Analysis of Gust-Response Characteristics of 
Vane-Airplane Combination 


In order to obtain theoretical verification of the apparent fact 
that the airplane response may be neglected in the calculation of gust 
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power іп the relatively short wavelength range, an analysis of the problem 
was made in terms of the theoretical frequency-response functions of the 
airplane. If the assumption is made that there are no applied controls 

or disturbances other than the gust, the following relationships for 

the rigid airplane can be written: 


7 n (8) 
ag & 


where the bracketed terms represent the frequency-response functions of 
the airplane which relate the angle of attack and rate of pitch to gust 
inputs. The performance of the vane as a gust-measuring instrument can 
be obtained by substituting equations (8) into equation (6) and solving 
for the ratio of the vane response to the gust input 


vU ee) ee) E > 


The airplane-frequency-response functions which appear in equa- 
tion (9) were calculated in an unpublished analysis of the two-degree- 
of-freedom response of the test airplane for the test conditions. The 
analysis included the effects of unsteady lift on the pitching moment 
and lift due to gusts as well as an approximation of the effect of a 
nonuniform distribution of angle of attack across the wing span due to 
short-wavelength gusts. With these airplane frequency-response functions, 
the frequency-response function of the vane-airplane combination as an 
instrument for measuring gust vertical velocity was calculated from 
equation (9) and plotted in figure h. From this figure it appears that, 
for the test configuration of vane and airplane, the vane performs very 
well as a gust-measuring instrument at frequencies above 1 cycle per 
second. Below l cycle per second it would seem that the vane measure- 
ments would be uséful, if corrected for the effect of the airplane 
response, down to approximately 0.1 cycle per second. Below this fre- 
quency the airplane response becomes the gréater part and vane measure- 
ments are of no further use for obtaining gust vertical velocity. 


Q 


The fractional error in the power spectral density of vane- ; 
indicated gust vertical velocity without corrections for the airplane 
response can be determined in terms of the performance function of 
equation (9) by the use of the theorem of equation (5) as follows: 
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зарби) enel... эш). [ur (20) 


PSD (wg ) PSD (we) x Oe 


where PSD( ) represents the power spectral density of the quantity in 
parentheses and wy, is the vane-indicated vertical velocity (Way = Vay). 


The error in the vane-measured power spectrum due to airplane response 
can be calculated from figure h and equation (10). The error is less 
than LO percent for frequencies greater than 1 cycle per second. 


PRESENTATION OF RESULTS 


The power spectral density of gust vertical velocity was calculated 
from measurements made with an angle-of-attack vane by following the 
procedure recommended in reference 4. Figure 5 represents the average 
power spectrum obtained from three test runs of from 50 to 60 seconds 
each made on two different days. The data approximate a straight line 
which has a slope of -2 when plotted on logarithmic paper. Each point 
has associated with it a band of confidence from the standpoint of sta- 
tistical reliability within which the true value lies for a given proba- 
bility (refs. 5 and 4). Тһе confidence band for 90-percent probability 
was calculated for the faired spectrum and is shown by the dashed lines 
in figure 5. Each calculated point represents the average power con- 
tained in the band width which is centered about it and extends halfway 
to the adjacent points. It should be noticed that the band width over 
which the average is taken is increased with increasing frequency with 
the result that the reliability is greater at the higher frequencies. 


The band widths over which the average power was calculated are 
illustrated better in figure 6. In this figure are compared the power 
spectra of gust vertical velocity calculated from an upwind test rum and 
a crosswind test run which were obtained on the same day only a few 
minutes apart. The results are plotted as bar graphs. The height of 
each bar or step represents the average power density for the band 
width over which it extends. The 90-percent confidence limits for fig- 
ure 6 are the same as those given in Pigure 7, which shows a comparison 
of the power spectra of gust vertical velocity and gust horizontal | 
velocity normal to the flight path recorded simultaneously. The 90-percent 
confidence limits, although shown related to the straight-line fairing of 
the vertical-velocity spectrum, are applicable to both. The assumption 
of no response of the airplane is made in the calculation of each 
spectrum. 


The standard deviations of gust vertical velocity for the upwind 
and crosswind test runs on the same day were 2.9 and 3.0 feet per second, 
respectively. The standard deviation for the vertical and horizontal 
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components of gust velocity recorded simultaneously were 2.1 and 4.8 feet 
per second, respectively. The comparatively high standard deviation of 


the horizontal velocity is a result of Low-frequency power in the record | 


which resulted from a lightly controlled Dutch roll oscillation of the 
airplane and is noticeable in the low-frequency values of figure 7. 


The probability distribution calculated for the 20-second trial 
record used in obtaining figure 5 is shown in figure 8 for comparison 
with the fitted normal distribution. 


DISCUSSION 


The power spectrum of gust vertical velocity is presented in fig- 
ure 5. For comparison, the lower frequency spectrum of reference 2 is 
also shown. Both spectra are compared with a straight line which has. 
a slope of -2. It is evident that atmospheric turbulence is well 
represented by a function which varies inversely with the square of the 
frequency. The spectrum form which has been found to be suitable for 
describing isotropic turbulence for measurements obtained in this manner 
is represented by the function (ref. 5) 


1 + 3(exfT)* 
(әкет, 


G(f) = а(0) (11) 


The constant T is the scale of turbulence which has the dimensions of 
seconds or feet depending upon whether f is frequency or inverse wave- 
length. The constant G(O) is proportional to both the scale of turbu- 
lence and the mean-square turbulence velocity. The high-frequency 


asymptote of equation (11) is proportional to f-2 which is in agree- 
ment with the experimental observations. The theoretical spectrum has 

a break frequency at the value for which the wavelength is equal to От. 
Below. this frequency the value of the power spectrum approaches the _ 
constant value G(O). None of the data of figure 5 extend to wavelengths 
long enough to confirm this characteristic. | 


The power spectrum of figure 5 is the average of three records, 
two of which were taken on the same day only a few minutes apart. One 
of the two records which were taken on the same day was taken while 
flying upwind; the other was taken while flying crosswind. Іп figure 6 
the power spectral densities of these two records are compared. The 
average wind direction appears to have no effect on the measurements. 


The spectra in figure 6 are presented as bar graphs in order to illustrate . 


better the band widths over which the average power is calculated. The 
same band widths are used in figures 5 to 7. 
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The third record was obtained on a different day having similar 
weather conditions. This record included simultaneous measurements of 
angle of attack and angle of sideslip. The power spectral densities 
of gust vertical velocity and gust horizontal velocity normal to the 
flight path are compared in figure T. The calculations for each com- 
ponent were made by assuming no airplane response. Тһе statistical 
confidence limits are shown for the faired vertical-velocity spectrum. 
Similar limits are applicable to the horizontal component. With allow- 
ances for the large low-frequency power in the horizontal-velocity 
spectrum which is the result of a lightly controlled Dutch roll oscilla- 
tion of the airplane, the agreement between the two components is good. 


Many problems concerned with gust disturbances involve considera- 
tions of normality, homogeneity, and isotropy of turbulent air. The 
probability distribution of gust vertical velocity was calculated only 
for the 20-second trial record which was used to determine the effect 
of neglecting the airplane response. Figure 8 shows the result com- 
pared with the fitted normal distribution. Figures 6 to 8 may encourage 
assumptions with regard to considerations of normality and isotropy. 
Figure 7 suggests that at least two-dimensional isotropy exists. The 
use of angle-of-sttack and angle-of-sideslip vanes in this fashion with 
the addition of a sensitive airspeed recording system for measuring | 
airspeed fluctuations due to gusts would be very useful for making 
detailed studies of isotropy of atmospheric turbulence. By measuring 
three components of gust velocity with respect to the airplane, theo- 
retical relationships between longitudinal and transverse correlations 
for isotropic turbulence could be checked. 


The results of this and other experimental determinations of the 
spectrum of atmospheric turbulence seem to justify the use of the 
spectrum form of equation (11) for those problems that require the 
determination of the effect of rough air on the performance of an 
aircraft control system. More long-wavelength data are needed in 
order to define the spectrum in the vicinity of the break frequency. 

The two primary characteristics of the turbulence spectrum are the scale 
and intensity. It must be determined over what range of values these 
characteristics vary and, from the standpoint of statistical probability, 
what design values should be chosen for specific problems. Іп particular, 
the value of the scale of turbulence is important in those problems for 
which the frequency range of interest is affected by the break frequency 
ОҒ the turbulence Spectrum. 


For the purpose of accumulating data on atmospheric turbulence over 
& wide range of wavelengths, flow-direction vanes of the type used in 
this investigation should be very useful except at the very low fre- 
quencies with respect to the response of airplanes. As was indicated 
in figure h and the section entitled "Method of Analysis, the angle-of- 
attack vane gives a direct measurement of the gust vertical velocity at 
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those frequencies for which the airplane response is negligible. At 
lower frequencies the vane indications can be corrected for the response 
of the airplane to afford a measurement of the gust vertical velocity. 
In this range of frequencies the gust power spectrum which is calculated 
from the corrected vane indications can be used to calculate the ampli- 
tude characteristics of the airplane frequency-response function or to 
check the derived frequency-response function. At still lower frequen- 
cies, the alrplane response quantities are the primary measurement; во 
the frequency-response function for the airplane-autopilot combination 
mist be used (since it seems desirable to stabilize the phugoid mode 

of the sirplane under these conditions). 


CONCLUDING REMARKS 


A flow-direction vane mounted on an airplane has been shown to 
provide a good measurement of gust velocity in the range of frequencies 
above the short-period response of the airplane. This use of vanes 
affords a simple, direct method of obtaining the power spectral density 
of atmospheric turbulence in that frequency range. By correcting the 
vane measurements for airplane response to gusts, the spectrum can be 
extended to Lower frequencies which include a large part of the short- 
period range. Therefore, when airplane response is measured for the 
purpose of calculating the power spectrum of gusts, an advantage is to 
be gained by recording angle~of-attack and angle-of-sideslip vane 
indications similtaneously. In addition to extending the range over 
which the spectrum can be calculated to higher frequencies, the spectrum 
calculated from the vane readings, when corrected for airplane response, 
can be used to check the results of the lower frequency measurements 
which require the use of theoretical airplane-response characteristics. 


The measured power spectral density of gusts displays an inverse- 
square relationship with frequency. That is, the average power per unit 
band width varies with the square of the wavelength over the frequency 
range covered. The wavelengths corresponding to the frequency range 
covered in these measurements are from 220 feet to 10 feet. This | 
inverse-square variation follows the trend exhibited in other measure- 
ments of gust spectra covering lower frequency bands. Until there are 
accumilated sufficient data on the power spectral density of atmospheric 
turbulence to set up design criteria based on statistical considerations, 
this form of the spectrum of turbulence should be useful for investigating 
the effects of turbulence over those frequency ranges of interest to 
airplane or control-system designers, provided the scale of turbulence 
can be assumed to be large. For those problems for which the frequency 
range of interest is affected by the break frequency of the turbulence 
spectrum, a value of scale mst be assumed for use with the theoretical 
spectrum. 
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The power spectral densities of gust vertical velocity and gust 
horizontal velocity perpendicular to the flight path were found to be, 
for practical purposes, equal. This result is suggestive of two- 
dimensional isotropy. Information concerning isotropy in atmospheric 
turbulence in the high-frequency range could be obtained readily by 
use of a sensitive airspeed recording system with angle-of-attack and 
angle-of-sideslip vanes for measuring three components of gust velocity 
with respect to the airplane. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., August 2h, 1954. 
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TABLE І.- SUMMARY OF OBSERVED METEOROLOGICAL AND TEST CONDITIONS 


Record length, вес. . e es ss.. 
Altitude, ҒҰҰ............. 
Free-sir temperature, °C ......- 
Wind direction . . . e. + s o «© «e s 
Approximate wind velocity, knots . . . 
Flight direction . . . a e a - š š 
Standard deviation of vanesiniicated 

gust vertical velocity, ft/sec ... 
Standard deviation of vane-indicated 

gust horizontal velocity, ft/sec . . 


va 
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Figure 1.- Installation of flow-direction vanes for measuring gust 
velocities. 
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Figure 2.- Symbols and axes used in analysis. Positive direction of 
angles and velocities is shown. 
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Figure 5.- Comparison of power spectral density of gus 
pre as obtained from corrected and uncorrected flow-direction vane. 


V = 220 feet per second. 


Phase angle ,deg 


Amplitude ratto 


Frequency ,срѕ 


Figure h.- Frequency-response function which describes the performance 


of yane-airplane combination as gust-measuring instrument, [v/a] - 


V = 220 feet per second. 
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Figure 5.- Average power spectrum of gust vertical velocity obtained 
from three records shown for comparison with data of reference 2. 
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Figure 6.- Comparison of power spectral density of gust vertical veloc- 
ity calculated from an upwind test run and a crosswind test run in 
the same area. V = 220 feet per second. 
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Figure 7.- Comparison of power spectral density of gust vertical veloc- 
ity and horizontal velocity normal to flight path recorded simulta- 
neously. V = 220 feet per second. 
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Figure 8.- Probability that ratio of deviation to standard deviation of 
gust vertical velocity will exceed a given value. 
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